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According to Bohmian dynamics, the particles of a quantum system move 
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^(x, t). We show that for simple one-dimensional systems any initial probabil- 

— ity distribution of a statistical ensemble approaches asymptotically |^(x,i)| 2 

■M | if the system is subject to a random noise of arbitrarily small intensity. 
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According to the de Broglie-Bohm quantum theory of motion [1-3], the n particles of a 
quantum system described by a wavefunction 

*(X, t) = R(X, t )e^ s ^\ X = {x u x 2 , ...Xi, x n } (1) 

that obeys the Schrodinger equation, follow trajectories whose velocity is given by the equa- 
tion 



1 



which can be solved knowing the initial conditions X(t — 0). Note that in this equation the 
actual value of the velocity (in configuration space) is the realization of a particular value of 
the velocity field 1/mV gS(X,t) at the actual position X(t) of the system. 

Consider now an ensemble of systems, associated with a probability distribution in con- 
figuration space, p(X,t). Both p and the quantum distribution probability, l^l 2 , obey a 
continuity equation [2], 

dg(X,t)/dt + V x -(g(X,t)X) = 0, <? = P,|*| 2 - (3) 

It is clear that if p and |\1/| 2 coincide at a given time t a , they coincide at all times (that is called 
the equivariance property) so that the predictions of quantum mechanics are reproduced. 
The condition p(X,t Q ) = \^(X,t )\ 2 is then sometimes taken as a postulate of the Bohmian 
quantum theory of motion (cf. e.g. [3]). However, this restriction of the theory to particular 
ensembles is not logically justified, since in the theory trajectories have an independent reality 
of their own [4] , and in fact there have been attempts to deduce the quantal distribution from 
the dynamics of the system. Precisely because of the existence of trajectories this problem 
can be considered to be akin to that of classical statistical mechanics [5]. Valentini [6,7] 
has deduced a "subquantal" H-theorem, which involves hypothesis similar to those adopted 
in the demonstration of the H-theorem in classical mechanics, assuming that the orbits 
of the system are sufficiently complicated, so that they sample all the accessible regions 
of configuration space. Indeed, many recent studies deal with the chaotic features of the 
deBroglie and Bohm trajectories for systems with several degrees of freedom (cf. e.g. [8]). 

We shall follow a different line of thought, inspired by the original paper by Bohm [9], who 
studied the consequences of random collisions among the electrons inside atoms or molecules, 
on the electron distribution. 



Let us suppose that the electrons of a given statistical sub-ensemble of molecules are 
initially described by the stationary wavefunction of the ground state, ip a (xi,t) , and by 
a distribution of particles, not necessarily corresponding to \ip a \ 2 - Let us further suppose 
that, during some amount of time, these molecules interact with a second set of uncontrolled 
molecules, described by the wavefunction (f>b(x 2 ,t). 

Initially (for t — > — oo) the wavefunction of the two colliding molecules is 

ty (x 1 ,x 2 , r 12 ; t) = ipaixi, t)(f) b (x 2 ; t)d (r 12 ; t), (4) 
while during the collision the wave function must be written as the linear combination 

^(x 1 ,x 2 ,r 12 ;t) = J]C a j3 )7 * aA7 (a;i;{a;2,ri2;t}) (5) 

with 

^ a ,/3, 7 (xi, x 2 , r 12 ; t) = ip a (xi; t)(j>/3(x 2 ; t)d 7 (ri 2 ; t) (6) 

where ip a (xi,t) represents a stationary state of molecule-1, ipp(x 2) t) represents a stationary 
state of molecule-2, and c? 7 (r 12 ,t) represents a state of relative motion of the center of mass 
of the colliding molecules (typically a wave-packet). The details of the initial state, that is 
^a(xi), (pb(x 2 ) and d (ri 2 ), determine, via the Schrodinger equation, the amplitudes C a ,p,y of 
the different stationary states of the compound wave function. 

Note that, in order to calculate the velocity of electron- 1, this wave function may be 
rewritten as the linear combination of eigenstates of molecule-1 

t) = J2 X«(x 2 (t), r 12 (t); t)Mxi; t), (7) 

a 

in which the coefficients X a (x 2 (t), ri 2 (t); t) incorporate the dependence in x 2 and ri 2 , and the 
summation over the indices /3 and 7. The velocity of the electron in the interior of molecule-1 
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is obtained using this wave function by means of Eq.(2), being thus sensitive to all the initial 
conditions, including, the positions of the electrons of molecule-2, (^(i — — °°)) an d the 
relative positions. Similar expressions may be adopted for the velocities associated to x 2 and 
ri2- 

After the collision (for t — > +00) the linear combination in Eq.(5) may be simplified, for 
all practical purposes, by retaining only one component, for example that of a = (3 = a in 
the case of elastic scattering. This is a kind of "natural" wave function collapse (see [2,10]), 
and guarantees that after the collision process the system is left in a stationary state. The 
question one would like to address is how the collision process affects the original distribution. 

In order to find an answer, we shall study the evolution of the particle distribution of one 
of the subensembles, that of coordinate aq, replacing the detailed description of the other 
variables (X2 and r^) during the collision process by the action of a random noise. The 
stochastic noise will represent the effect of the deviation of all those other degrees of freedom 
from their average value. 

We shall solve the equations of motion, studying numerically the resulting distribution 
function of the electrons, instead of making the strong statistical assumptions adopted by 
Bohm [J Our aim will be to show that the system tends to a distribution described by \^\ 2 , 
starting from an arbitrary initial distribution. 

We first study some properties of the time evolution of the distribution of particles 
governed by Bohmian dynamics, in relation with its response to a stochastic perturbation. 

1 Bohm's proof was based on the assumption of a continuity equation (Eq.(3) in [9]), even if the 
random perturbation acting on the electrons inside the molecules would require the more general 
Fokker-Planck equation, including sources and diffusion terms. 
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In particular we consider the effects of stochastic noise on one-dimensional, non-degenerate 
systems in non-stationary states. Starting from arbitrary initial ensemble distributions, we 
will show that they tend to the quantal distribution. 

In a realistic calculation, the noise intensity should be determined in accordance to the 
coefficients of the linear combination (7), and as a consequence should be dependent on the 
position along the trajectory. It would be natural to expect that during the reaction the 
noise will reach a peak value around the distance of closest approach between the colliding 
molecules, tending to zero for t — > ±00. In the following, we shall limit ourselves to a very 
simplified description, solving the equations of motion (the one-dimensional equivalent of 
Eq.(2)) 

x i = \/mVS{x i ,t) + Cri i , i = 1, 2, ...TV (8) 

where % labels the N elements of the statistical ensemble of one-particle systems, initially 
distributed according to a given function po(x) = p(x,t = 0), and subject to a random 
white noise rji(t), whose intensity is controlled by the parameter C (in what follows C will 
be measured in hundredths of the speed of light, which is a natural scale for the studied 
systems). The ^-independent character of the noise that we use, would produce a complete 
diffusion in the whole space. This is avoided by confining the system in a box with infinite, 
perfectly reflecting walls. 

Let us thus numerically study the case of a statistical ensemble of N = 5000 test systems, 
each of them representing an electron confined in a box of length L = 2A with wave-function 

*(x, t) = a 1 (j> 1 (x)e~ iElt/h + a 2 Mx)e~ tE2t/ \ (9) 

where <f>i and 02 are the two first stationary states of the well, taken with amplitudes a±, a 2 , 
which in a more general case would represent typical values during the reaction process (cf. 
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Eq.(7)). Here for simplicity we shall use the constant values a\ = a 2 = l/y/2. We have 
solved the stochastic differential equation (8) using Heun algorithm [11]. 

We shall compare the results obtained with two different initial distributions, a uniform 
one, po(x) = p un (x) = 1/L, and another coinciding, at t — 0, with the usual quantum 
distribution, po(x) = p q {x) = \^(x, 0)| 2 . The time evolution of po(x,t) starting from p un (x) 
is shown in Fig.l for the case of C = 0.1, which is a rather large noise intensity, corresponding 
to about 2% of the average velocity . It is seen that the approach to \^(x, t)\ 2 is already 
significant after one oscillation of the distribution. 

We shall measure the difference between the distribution p(x, t) evolving from the initial 
one, po(x), and the quantal distribution \^(x, t)\ 2 making use of the function 

X%(t) = J \p(x,t)-\*(x,t)\ 2 \dx. (10) 

In Fig. 2 we present the results for three different values of C. For each one we present 
the time evolution of x 2 ■ To be noted that for both initial conditions (p un (x) and p q (x)) the 
systems evolve under the influence of the random noise; therefore, even starting from p q , x 2 Pq 
is in general different from zero. Similar results are obtained using other measures of the 
difference between p(x,t) and ^/(x, t), in particular the relative entropy used in [7,8]. 

It can be seen that y 2 decreases until it reaches a minimum value, around which it 
stabilizes. This value coincides (apart from statistical fluctuations arising from the finite 
number of particles) with the one attained by x 2 q ■ This feature is seen in all the three 
cases presented, although the time needed for x 2 to stabilize is longer (r= 60, 600 and 6000 
periods) for smaller values of C (C— 0.01,0.001 and 0.0001). It is also important to note that 
the asymptotic value of x 2 becomes smaller for smaller C, which means that the obtained 
particle distribution is closer to the quantum mechanical value \^(x, t)\ 2 the smaller the value 
of C is. This is illustrated by the distributions presented on the right side of Fig. 2. It can 
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be seen that reducing the noise, x 2 tends to an asymptotic value that is not equal to zero, 
what is mainly (90%) due to the finite statistical sampling (5000 systems distributed over 
200 x-subintervals), and in a much smaller extent to the unavoidable errors in the numerical 
integration of the equations of motion (performed with At = 1/500 of a period). In fact, 
we have verified that eliminating the noise, that is, putting C equal to zero, the asymptotic 
value is equal to 0.20, very close to the value obtained for C = 0.0001. 

In order to obtain some insight in the results discussed above, it is convenient to rewrite 
the Bohmian equation of motion as a second-order newtonian equation for the acceleration [2], 

mxi = -VU d - VU qu , (11) 

where U c \ is the classical potential and 

U qu = -(k 2 /2m)V 2 R/R (12) 

is the so-called quantum potential. 

In the case of the particle in the box U d = and thus, aside from the stochastic noise, 
the trajectories are determined by the behavior of U qu . This is a time dependent potential, 
which goes to infinity where R = 0, that is at the nodal points. Instead at the points where 
R has maxima, U qu is almost flat. 

In the absence of noise, the Bohmian trajectories may be thought as those of classical 
particles moving in the time dependent potential U qu , with the condition that the initial 
velocities, needed for integrating Eq.(ll), must be those determined by the wave-function of 
the system (cf. Eq.(2)). 

In our case, with only two eigenstates in the linear combination of ^ (cf. Eq.(9)), the 
density oscillates in time with Bohr's period T = 2nh/(E 2 — Ei), and thus the potential 
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U qu (x,t) and the trajectories are also periodic (see Fig. 3 and Figs. 4a,b) Q. The trajectories 
near to the nodal points are compressed and re-expanded periodically by the strongly varying 
quantum potential. Such a compression-expansion is instead very soft for the trajectories 
close to maxima of R (note R in the denominator in Eq.(12)). 

The effect of the stochastic noise is to break this periodic pattern (see Fig. 4c). This 
breaking can be very effective, because any slight perturbation in the coordinate during the 
compressed phase makes the particle dramatically change its subsequent trajectory in the re- 
expansion process. The net effect of all these changes of trajectory is a higher residence time 
around the maxima of R, that produces a distribution of particles resembling |\I/| 2 in spite 
of the initial distribution. We then see why the distribution |\I/| 2 displays the remarkable 
stability against random perturbations of the system found above (cf. Fig.2). 

We have checked that the above results are neither a special feature of the chosen po- 
tential, nor of the particular considered wave-function, by repeating the calculation with 
potentials of different shapes inside the walls and with different linear combinations of eigen- 
states, obtaining always a similar behaviour for the density distribution. In the case of the 
parabolic potential we have also analyzed the classically forbidden region, that is the region 
where the classical potential is larger than the average energy of the considered wave func- 
tion, and found again that the density distribution approaches |\l/| 2 as the noise intensity is 
set to smaller values. Let us however make some specific comments on the case of stationary 
states. 

According to Bohmian dynamics, in the case of a stationary state (^(x, t) = <p(x)e~ tEt l h ) 

2 In general, for more complex linear combination of states of the box, the periodicity will be 
controlled by the ground state frequency, of which all the other frequencies are multiples. 
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the electrons are at rest (cf. Eq.(2)), and the addition of a stochastic noise would give 
rise to a pure diffusion process, leading inevitably to a uniform distribution inside the box. 
However this situation is not physically meaningful since, as stated above, the stochastic 
noise represents the (random) deviations of the single-particle like wave-function in eq.(7) 
from its average linear combination, which for a pure stationary state is just zero. In other 
words, the non-stationary character of the wave-function and the stochastic noise go together. 

We conclude that under the action of random noise, at least for very simple, one- 
dimensional systems, an ensemble of particles, initially arbitrarily distributed, and governed 
by Bohmian dynamics progressively loses the "memory" of the initial distribution as time 
proceeds, and approaches the distribution \^(x, t)\ 2 , a process which is based on the different 
stability of the underlying Bohmian trajectories. The origin of the random noise may be 
naturally found in the uncontrollable character of the position of the particles in Bohm's 
theory. 
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Figure captions 



Fig. 1 

Evolution of the particles distribution during the first oscillation corresponding to an 
initially uniform distribution subject to a noise intensity of C — 0.1 (continuous line) and 
the corresponding |\l/(a;)| 2 (dashed line). 

Fig. 2 

On the left hand side the evolution of Xp un (t) (solid line) and xj> q (t) (dashed line) is 
shown as a function of time for three different values of C (0.01, 0.001 and 0.0001). On the 
right hand side the distribution p(x, r) evolving from p un (solid line) for the three different 
C values at a fixed values of time (r= 60, 600 and 6,000 periods, respectively), is compared 
to \^(x, t)\ 2 (dashed line). 

Fig. 3 

The periodic, time-dependent quantum potential U qu (x,t) (in eV) is shown as a function 
of x (in A) and t (in periods). Note the nodal point appearing as peaks in U qu (x,t). 

Fig. 4 

a) Trajectories of a set of particles initially distributed according to \^(x,t = 0)| 2 , and 
without any noise, b) Trajectories of a set of particles initially uniformly distributed (p G = 
p un = 1/L) without any noise, c) The same as in b) but under the action of the stochastic 
noise (C = 0.1 in this example). 
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